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Abstract 

In this paper we present a construction of a new class of explicit solutions to the 
WDVV (or associativity) equations. Our construction is based on a relationship be- 
tween the WDVV equations and Whitham (or modulation) equations. Whitham equa- 
tions appear in the perturbation theory of exact algebro-geometric solutions of soliton 
equations and are defined on the moduli space of algebraic curves with some extra 
algebro-geometric data. It was first observed by Krichever that for curves of genus 
zero the r-function of a "universal" Whitham hierarchy gives a solution to the WDVV 
equations. This construction was later extended by Dubrovin and Krichever to alge- 
braic curves of higher genus. Such extension depends on the choice of a normalization 
for the corresponding Whitham differentials. Traditionally only complex normalization 
(or the normalization w.r.t. a-cycles) was considered. In this paper we generalize the 
above construction to the real-normalized case. 



Introduction 

In the beginning of the 90's, while studying deformations of 2D topological field theories 
(TFT), E. Witten, R. Dijkgraaf, E. Verlinde, and H. Verlinde ( ||DVV91b| , PVV914 |Wit9q| ) 
wrote down the following overdetermined system of non-linear PDEs for a function J-{t) = 
J-(to, t±, . . . ): 



(WDVV) ^ja^oAm) -1 ^ = ^a(^oam)- 1 ^ 



where T a = -= — , and the matrix r] a p = T§ a p is constant and non-degenerate. These equa- 



tions are now called the WDVV equations. They appear in the following way ||Dij98 



One can show that the structure of a 2D TFT is equivalent to a structure of a Frobenius al- 
gebra A, i. e., a commutative associative algebra with a unit and a symmetric non-degenerate 



bilinear form r\ such that r)(a* b, d) = r/(a, b*d). Let {4> a } be a basis of A (<f> a correspond to 
primary fields in TFT) and let j] afi = r}((f) a , <f>p) = (fefo?)o> ^7 = '/(".,- o i*Q- ) = (0 a 0/3</> 7 ) o - 
Deformations of a TFT structure considered in [ DV V91b , DV V91aj , |Wit90|| correspond to 



potential deformations of a Frobenius algebra, i. e., there should exist a function .F(t), called 
the WDVV potential, such that 770-/3 = doapj 7 ^) and c a/ g 7 (t) = d a/ g 7 .F(£) . Then the WDVV 
equations are just the associativity conditions for the deformed algebra structure. 

It is now clear that the theory of the WDVV equations and its differential-geometric coun- 
terpart, the theory of Frobenius manifolds, are related to a whole spectrum of applications 
ranging from classical differential geometry (Darboux-Egoroff metrics, n-orthogonal curvi- 
linear coordinate systems, deformations of singularities) to quantum cohomology, Gromov- 
Witten invariants, integrable hierarchies, and the Seiberg-Witten equations. 

Solutions of the WDVV equations that can be obtained from the theory of Whitham 
hierarchies correspond to the topological Landau- Ginzburg theories and minimal models. 
For A n minimal models, a Frobenius algebra structure is defined on the space of degree n — 2 

polynomials in p with the help of the superpotential W(p) = — by 

n 

u * v = u(p)v(p) mod W'(p), 

1 v u(p)v(p) 
(uv) = -re Soo ^^-dp. 

A Frobenius algebra structure is then deformed by deforming the superpotential W, 

W(p\a) = h a n _ 2 p n ~ 2 + h dip + a , 



where we used the notation W(p|a) to separate the variable p from the deformation param- 
eters a n -2, ... , ao- Let t(p) be an n th root of W(p|o) in the neighborhood of infinity, 

r (p) = nW(p) = p n + • • • , t{p) =p + 0{p- 1 ). 

Then the deformed basis of primary fields <p a is given by 

a) ^=-L%»^ t .., 

a + 1 dp 

where Q a (p) = [t a {p)] is a principal (i.e., polynomial) part of the Laurent expansion of 
t a in the neighborhood of infinity. The flat coordinates t a on the space of the deformation 
parameters are given by t a = — resoo t~( a+1 ^pdW, and (j) a {p\t) = —d ta W(p\t). Then 

Va/3 = (0a, 4>I3) = ~ ^oo ^/A, = 0~a+f3,n-2, 

> dVt a+ idVLp + idVl 1+ i 
c a ^t) = (coeff) ^ res, — . 

q s \dW(q s )=0 

In ||DVV91b| it was shown that c aj a 7 (t) satisfy the integrability conditions and therefore 
c a (3~f(t) = da^J-'it). The closed expression for this WDW potential jF(rj) was identified by 
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I. Krichever ||Kri92|| with the logarithm of the r-function of a certain reduction of the genus 



zero Whitham hierarchy. 

Whitham equations (or modulation equations) appear in the theory of perturbations 
of exact algebro-geometric solutions of soliton equations and describe "slow drift" on the 
moduli space of algebro-geometric data. These equations can be defined with the help of 
certain differentials dQi on the universal curve, each differential generating a corresponding 
Whitham flow on the moduli. A large class of solutions of the Whitham hierarchy is given 
by the so-called algebraic orbits. Such solutions depend only on finitely many parameters 
and can be constructed as follows. One starts with a finite-dimensional moduli space, called 
the universal configuration space, that consists of a curve T, punctures P a , and a pair of 
Abelian integrals E and Q. Then one defines special coordinates on this space (Whitham 
times), picks a leaf defined by fixing some of them, and maps this leaf to the moduli space 
of algebro-geometric data in such a way that coordinate lines go to Whitham flows. All 
information about such algebraic orbit can be encoded in a single function r(t) that depends 
only on the moduli. This function is called the r-function of an algebraic orbit. 

In the genus zero case Abelian integrals become polynomials, and if we choose an algebraic 
orbit with Q = p, then E can be identified (up to normalization) with the superpotential 
W, Qi define a basis of the corresponding Frobenius algebra, Whitham times tj give flat 
coordinates on the orbit, and E(t) = lnr(t) is a WDVV potential. 

This approach can be generalized to moduli spaces of curves of higher genus. One new 
feature of a higher genus case is a more complicated topology of the moduli space. As a 
result, the hierarchy has to be extended to include certain (multivalued) differentials dVtj^ 
that generate additional Whitham flows. Another important issue is a choice of the nor- 
malization. Namely, in the genus zero case, differentials dQi are completely determined by 
their expansions in the neighborhoods of marked points P a . In the higher genus case, these 
conditions define dQi only up to a holomorphic differential. This ambiguity is fixed by intro- 
ducing a normalization condition. There are two main choices — real normalization, which 
is defined by 
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dQi 



o Vce#i(r s ,z) 



and complex normalization (or normalization w.r.t. a-cycles). Complex normalization re- 
quires making a choice of a canonical basis B of cycles in the homology of T g and is then 
defined by 



dQi = Va fc G B. 



Whitham equations were originally derived in [|Kri8§|| for real-normalized differentials. How- 



ever, after the relationship between Whitham equations and WDVV equations was found 
Kri92|| and |Pub92|| , the focus shifted to the complex normalization condition |[Kri94|| . 



in 



There were two main reasons for this. First, in the complex-normalized case the r-function 
is holomorphic, which is important for string theory applications. Second, the derivation of 
the expression for the r-function relied on the Riemann bilinear identities. Corresponding 
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identities for real normalized differentials are technically more complicated. However, com- 
plex normalization has certain disadvantages. In particular, it is well-defined only on the 
extended moduli space that incorporates the choice of a canonical basis B into moduli data. 

In this paper we develop a real-normalized version of the above approach. In this case, 
Whitham hierarchy can be defined on the usual moduli space of curves with some extra 
algebro-geometric data. First, we define a real leaf in the universal configuration space, 
introduce Whitham coordinates on this leaf and map it into the moduli space in such a 
way that the resulting differentials are real-normalized. Then we prove the real-normalized 
version of the Riemann bilinear identities (generalized to multivalued differentials). Using 
these identities we find the formula for the r-function of an algebraic orbit and prove the 
following theorem for F(T) = lnr(T). 

Theorem. The third derivatives of F(T) are given by the following formula 



d TATBTc F{T) = 3? 



dflAd^Bdflc 



^ TeSqs dEdQ 

q a \dE{q a )=0 ^ 



This theorem then implies that if we consider a special class of algebraic orbits by choos- 
ing dQ to be a real- normalized differential with a pole of order two at a puncture Pi, then 
Whitham flows define potential deformations of a Frobenius algebra structure with observ- 

ables corresponding to — — , and we have the following theorem. 

dQ 

Theorem. The logarithm of a r-function of a (reduced) real-normalized genus g Whitham 
hierarchy is a solution to the WDVV equation. 

Acknowledgments The author is grateful to I. Krichever for suggesting this problem and 
for constant attention to this work. The author also thanks Yu. Volvovski for many useful 
discussions. 



Whitham equations 

Whitham equations first appeared in the theory of perturbations of exact algebro-geometric 
solutions of soliton equations. Such solutions are defined by linear flows on the Jacobian 
Jac(r 9 ) of an auxiliary algebraic curve T g and are expressed in terms of theta functions. 
Perturbing these solutions by the so-called non-linear WKB (or Whitham averaging [|Whi74j| ) 
method ( ||FFM80| , |DM82| , |DN83|| ) results in a "slow drift" on the moduli M of algebro- 
geometric data. Equations describing this drift are called Whitham equations. 

Although Whitham equations are equations on the moduli, they can be conveniently 
written with the help of certain Abelian differentials dfi^(P|J) 



P G T g defined on the 



universal curve Mg, 



M. 
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Each of the differentials cLVLa is coupled with a corresponding Whitham time Ta defining 
A th Whitham flow on M.. Then, after making a special choice of connection on Mg, Whitham 
equations can be written in the following implicit form, 

8a(1Qb = dsd^A- 

This form of Whitham equations was first observed by Flashka, Forest, and McLaughlin 



TM80] for the KdV equation and hyperelliptic spectral curves. It was later justified by 



Krichever [|Kri88|l in a more general setting of (2 + 1) equations and general spectral curves. 



Whitham hierarchies 

It turns out that one can construct a whole hierarchy of Whitham equations. Following 
Krichever, we use algebro-geometric approach to construct Whitham hierarchies (Hamilto- 
nian approach to the theory of Whitham equations was developed in ||DN83|| ) . 



We begin with a local definition. Let T = {Ta}a&a De a collection of (real or complex) 
times, indexed by some set A, let z G D C C, and let {Qa(z\T)}a£A be a collection 
of functions, meromorphic in z, each Qa is coupled with the corresponding time Ta- The 
functions VIa(z\T) should be thought of as pull-backs via Whitham times of Abelian integrals 
VLa = J d^A from a leaf spanned by Whitham flows in A4, with z corresponding to a local 
coordinate along the fiber. Define a 1-form u by 

uj = J2^a(z\T) dT A . 

A 

Then, by definition, Whitham hierarchy is given by the generating equation 

(2) 5u A 5lu = 0, 

where 

Slu = V d B VL A dT B A dT A + ^ dz A dT A . 

, dz 

A 



This is equivalent to the following two equations, 

(3) £ s^ B ' c > (^f) (d B n c ) = o, £ e^ D Hd A n B )(d c n D ) = o, 

{A,B,C} ^ ' {A,B,C,D} 

where we sum over all possible permutations of a fixed collection of indexes and e is a sign 
of a permutation. 

Usually we have one marked index Aq e A with 

X = T Ao , p(z\T)=Q Ao (z\T). 
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Then, as long as — 7^ 0, we can make a change of coordinates from (z, T) to (p, T). In these 
coordinates, equation @ written for A , A, B takes the form 



(4) 



d A Q B (p\T) - d B tt A (p\T) + {Q A , Q B } (p\T) = 0, 



where 



{n A ,n B }( P \T) = d x n A 



dVt B 



- d x ^B 



dVt A 



dp 



dp 



is the usual Poisson bracket. Equations in this form are called zero- curvature equations. 
They can be interpreted as a compatibility conditions for the following Hamiltonian system, 



Then, if zero-curvature equations are satisfied, there exists (locally) a solution E = E(p\T) 

dE 

of this system and, as long as — — 7^ 0, we can again change coordinates from (p,T) to 

dp 

(E,T). In these coordinates, the above system takes the form 



Therefore, there exists (locally) a function S = S(E\T) such that Q A (E\T) = d A S(E\T). 
This function S(E\T) is called a prepotential of the Whitham hierarchy. 

Algebraic orbits 

One can use the above formalism to construct certain exact solutions of Whitham equations. 
These solutions depend on finitely many parameters and are obtained as follows ||Kri94|| . First 
we construct Whitham flows on finite-dimensional submanifolds, called algebraic orbits, of 
M. (note that Ai is usually infinite-dimensional) together with Whitham equations that they 
satisfy. Then these equations are extended to the whole Ai. By making a correct choice of 
algebraic data, one can obtain usual Whitham equations of the soliton theory. 

We illustrate this idea in the special case of the so-called dispersionless Lax equations. In 
this case we take algebraic curve of genus zero, T ~ CP 1 , with a single puncture P\ = 00. 
The moduli space is equivalent to a moduli space of local coordinates around Pi, 



d A E(p\T) = {E,Q A }(p\T). 



d A p(E\T) = d x n A (E\T) 



and its compatibility conditions can be written as 



d A n B (E\T) = d B n A (E\T). 



M=M 



0,1 — 



{r ~ CP 1 ; P x = 00; z{P)} ~ {u a , s 



I,---}, 



where 



z 1 (p) = t(p) = p + vip 1 + v 2 p 2 H , 
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and p is the standard coordinate on C. Note that in this case f2i = p, i.e., the marked 
index is i = 1. Then, by definition, the dispersionless Lax hierarchy (or dispersionless KP 
hierarchy) is a set of evolution equations on v s = u s (T): 

d^{p\T) = {t, (f ) + } = {t,fii} , where d t = 

This hierarchy can be thought of as a quasi-classical limit of the usual KP hierarchy, and al- 
gebraic orbits correspond to n th order reductions of the KP hierarchy, i.e., nKdV hierarchies. 
Namely, let 

E(p) =p n + u n . 2 p n ~ 2 + --- + u 

and define t(p) by the condition t{p) n = E{p). Then u s = p s (uq, . . . , u n _2) and we obtain 
a finite-dimensional leaf in M.. Corresponding evolution equations on U{ = Ui(T) can be 
written as 

diE{p\T) = + | in (p,T) coordinates, 

dip{E\T) = d x tti{E\T) in (E, T) coordinates. 

Then the solution E = E{p\T) (i.e., Ui = Ui(T)) can be obtained as follows ||Kri88| , [Kri92 . 
Theorem. Let 

(5) S(p\T) = Y,T&i(l>\u) = J2 T ^ + °( r1 )- 

i i 

Require that 

dS dE 

(6) l~(Qs) = for all q s such that —^(q s ) = 0. 

dp dp 

Equation (|6|) is equivalent to a collection of equations F k (u,T) = that implicitly define 
Ui = Ui(T). 

Then u(T) is a solution of the dispersionless Lax equations. 

The theorem follows from the fact that S defined above is a global prepotential, i.e., 
d t S(E\T)=Q t (E\T) 



r + 0(l _1 ) near P 1 
has no other poles 



From the definition of S it is clear that the first condition is satisfied. The role of the 
condition (^) is to ensure that diS(E\T) is holomorphic on T — Pi. 

Note that the condition @ can be also written in the form dS = QdE for some polynomial 
Q(p\T). Conversely, we can recover Tj from S and t by the formula 

Tj = — - resp 1 1 dS. 
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This observation motivates the following alternative approach. Consider the so-called uni- 
versal configuration space 

M (n, m) = {T = CP 1 ; P 1 = oo; [z] n ; E, Q) ~ {u , u n _ 2 ; bo,..., b m } , 

where 

E(p) =p n + u n . 2 p n ~ 2 + --- + u , Q(p) = b mP m + ■ ■ ■ + b , 

[z] n is an n-jet of a local coordinate near P x and we require that E = z~ n + O(z) near P x . 
By definition, put 

(7) dS = Q dE, Ti = ~ res Pl z l dS. 

Then 

didS(E\T) = dQi(E\T), 

condition fl7|) is equivalent to a collection of equations Tk = Tk(u,b), which can be inverted, 
Ui = Ui(T), bj = bj(T), and therefore we obtain a map from M.0,1 to an algebraic orbit in 

•Mo- 
T-functions and solutions to the WDVV equations 

For each algebraic orbit constructed above corresponds a so-called r-function. By definition, 

]nr(T) = F(T), 

where 



and dS is a (differential of) the prepotential of the corresponding algebraic orbit. Then we 
have the following theorem (Krichever |[Kri92| ) . 

Theorem. The derivatives of F(T) are given by the following formulas: 

diF(T) = re Soo f dS 
dijF(T) = reSoo t J dQi = reSoo t J dQi 

sr^ dVtidVtjdVt k dS 

d ijk F(T)= ^ res, dQdE ' where Q=dE- 

q s \dE(q s )=0 ^ 

Corollary. If we choose 

Q{p\T) = p, E(p\T u T n+1 ) =p n + u n _ 2 (T)p n - 2 + ■■■ + u (T), 
then F(T) is a WDVV potential for A n _\ Landau- Ginzburg model defined by a superpotential 

W(p\t ,...,t n - 2 ) = -E(p\t , t ±,...,—,0, ' 



n I n n + 1 
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Higher genus case 

The above approach can be generalized to algebraic curves of higher genus. This question 
was considered by Dubrovin |[Dub92| for Hurwitz spaces and by Krichever [|Kri94|1 in general. 
In the higher genus case, the "universal" moduli space 



[P)} 



M g , N = {r g ;P a ,r a \p) 

consists of the following algebro-geometric data: 

• smooth algebraic curve T g of genus g, 

• collection of marked points P a G T g , a = 1, . . . , N, 

• local coordinates z a (P) in the neighborhoods of P a , z a (P a ) = 0. 

For simplicity, we concentrate exclusively on a single puncture case, N = 1. 

To construct a realization of a Whitham hierarchy on M. g ,i, one can proceed as follows 



sec 



[Kri94|| ). Instead of polynomials f2j one has to consider Abelian differentials of the 
second kind dQi on T g with prescribed behavior near the puncture, 



(9) 



dtii = d(r + 0(t -1 )) near Fx. 



However, condition @ specifies dVti only up to a holomorphic differential. To define dfli 
uniquely, one has to impose certain normalization conditions. There are two choices — real 
normalization and complex normalization (or normalization w.r.t. a-cycles). Real normal- 
ization is defined by the condition 



(10) 
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o Vcg #i(r 3 ,z). 

To define complex normalization, one has to first choose a canonical homology basis for T g , 

B = {ai, . . . , a g ; b l , . . . , b g \ a h bj G #i(r 9 , Z), a { ■ aj = b { ■ b i = 0, a { ■ bj = 5 i:j } . 
As a result, it is necessary to consider the extended moduli space 

A; tl = {T g] p 1] z{p) ] B} i 

which is a covering of M. g ,i- The differentials dfli are then uniquely normalized by the 
condition 



dQi = Vat G B. 



As before, the marked index is i = 1 and p(P) = J P d£li 
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Universal configuration space 

In order to construct solutions corresponding to algebraic orbits, it is convenient to use 



the universal configuration space approach. Following [|KP97| , |KP98|] , define the universal 
configuration space for Whitham hierarchies, 

M g (n, m) = {T g ; P x ; [z] n ; E, Q} , 

to be the moduli space of the following data: 

• genus g Riemann surface T g , 

• marked point Pi G T g , 

• an n-jet [z] n of local coordinates z(P) in the neighborhoods of Pi, z(P\) = 0, 

• an Abelian integral E with pole of order n at Pi such that 

(12) dE ~ d(z~ n + O(z)) near Pi, z G [z] n , 

• an Abelian integral Q with pole of order mat Pi. 

More precisely, by Abelian integrals we mean pairs E = (dE, P<f ), Q = (dQ, P Q ), where dE, 
dQ are meromorphic differentials of the second kind on T g , holomorphic on T g — Pi, and 
with poles of order n+1 and m + 1 at Pi, and 



E(P) = f dE, Q(P) = [ 
Jpz Jp 



p 

dQ. 



Alternatively, one can choose a local coordinate z G [z] n and define E and Q to be pairs 
E = (dE, ce), Q = (dQ, cq), where near Pi, 

E ~ z~ n + c E + O(z), Q ~ c Q m z- m + ■■■ + cfz~ l + c Q + O(z). 

Note that after E is chosen, there is a preferred local coordinate z* G [z] n defined by 

E = z~ n near P x . 

The moduli space M. g (n, m) is a complex manifold of dimension 5g + n + m with at most 
orbifold singularities. We also need to consider smaller (4g + n — l)-dimensional moduli space 

M g (n) = {T g ;P 1 ;[z] n ,E}, 

as well as the corresponding extended moduli spaces J\A*(n,m) and Ji4*(n). 
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Construction of algebraic orbits 

A general approach for constructing algebraic orbits in the higher genus case is the following. 
First, choose a leaf V(n, m) in M* g (n, m). The precise definition of such leaf depends on the 
choice of the normalization. Then introduce special coordinates Ta, called Whitham times, 
on V(n, m), and define the prepotential S of the leaf by the formula 

dS = QdE. 

Whitham differentials dflA are then obtained from dS by 

d TA dS(E\T) = dn A (E\T), 

where we use a connection on Mg given by choosing E = const to be horizontal sections. 
Note that dS = J2 A T A dQA- The leaf V(n, m) is then mapped to the corresponding algebraic 
orbit 0(n) by the following sequence of maps: 



(13) 



V(n,m) ► V{n) ► 0{n) 

U U 

V°(n) 0°(n) 

{Tg-Py^zUE^Q} > {r g ;P 1 ;[z] n ;E} > {F g ; P ± ; z*} , 



where V°(n) is a leaf of a subfibration of V(n) defined by fixing the periods of dE. Then, in 
the coordinates (E,T) on the "universal curve" J\f g i over 0°(n), the differential dS(E\T) = 
Q(E\T)dE satisfies the equations dAdS(E\T) = dVt A {E\T). Therefore, dS is a prepoten- 
tial, 8 x Qa(E\T) = 8t a p(E\T), and we obtain solutions to Whitham equations. Note that 
Whitham coordinates on V(n, m) go to Whitham flows on the algebraic orbit. 

For each Whitham derivative &t a we define a dual "integral" operator § T in such a way 
that 

(14) I d TB dS= I dtt B = I dtt A = I d TA dS VT A ,T B . 

The exact definition of such operators § T can be rather non-trivial and usually includes 
integration with some weight over a cycle and, maybe, some correction terms. The above 
identities can be thought of as a generalization of Riemann bilinear identities for Whitham 
differentials, and they are proved along the same lines. The r-function of an algebraic orbit 
is then defined by r(T) = e F<yT \ and 



(15) F(T) = lj2 T * ( f dS - 

It encodes all information about the differentials dVL A and the curve T g , and its first and 
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second derivatives are given by 

i r 1 ^ r 



~ <f dS+\y^T B I dQ B = <f dS, 

1 Jt a 1 Tb Jt a Jt a 

d TATB F(T) = I dQ B = I 

JTa JT r 



dVtA- 

>T A JT B 

Note that the equality of mixed partials of F(T) corresponds to equation (fLip . Third deriva- 
tives of F(T) are then given by the residue type formula. 

Complex-normalized case 

In this case, one chooses 

t k = —-r resp, z k Q dE, k = 1, . . . , n + m, 
k 



th,k = <P QdE, t E)ak = (h dE, t Q: a k = (h dQ, t E ,b k = f dE, t Q , bk = (h dQ, 

J a k J a k Ja k Jb k Jb k 

as coordinates on Ai g (n,m). The notation a k indicates that the integral has to be taken 
over the right side of the cycle. A complex (3g + n + m)-dimensional Vc{n, m) C A4*(n, m) is 
defined by imposing the a-normalization condition on the Abelian differentials dE and dQ, 

dE= I dQ = \/a k G B. 

>a k J a k 

Note that the projection of Vc{n, m) to M. g (n, m) is not well-defined. Whitham coordinates 
on the leaf are then given by 

T k = — — resp. z k dS, k = 1, . . . , n + m, 

k 



Th,k = f dS, T E)k = ~ f dQ, T Qtk = <b dE, k = l,...,g 

J a k Jb k Jb k 

Corresponding Whitham differentials are divided into the following four groups. 

• dT k dS(E\T) = dQ k (E\T) are the usual complex-normalized Whitham differentials with 
prescribed behavior near the puncture, 

dQ k = d(t k + C^r 1 )) near P u I dQ k = 0. 

Jai 

• dr h k dS(E\T) = dQh,k{E\T) are holomorphic differentials that are dual to the canonical 
basis of cycles, 



/ 

J a, 



dSlhk — Sik- 
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• 8t e k dS(E\T) = dQ E ,k{E\T) are a-normalized and holomorphic everywhere on T g ex- 
cept for the a-cycles, where they have jumps, 

dtt E)k {P+ - P-J = S kl dE(P ai ), I dn Etk = 0. 

• Otq k dS(E\T) = dQQ } k{E\T) are a-normalized and holomorphic everywhere on T g ex- 
cept for the a-cycles, where they have jumps, 

dn Q , k (P+ - P-J = 5 kl dQ(P ai ), I dQ Qyk = 0. 
The duals § T are then given by 



<p dQ = res Pl Vdfl, <p dQ = — 



h.k 



dQ 



dn = — 

T E , k 27T. 



-1 L 



EdQ 



where the notation , dQ indicates that a certain correction terms have to be added to 
make the integral independent on small cycle deformations. In particular, 



(f dS= (f dS + T Etk E(atnb k ). 



This remark becomes very important in the real-normalized case and is discussed at length 
in the next section. 

The third derivatives of F(T) are given by the following theorem ( | |Kri94|| ) . 



Theorem. 



& 



T A T B T C 



F{T) = 



res. 



dVLAdVLsdVl 



c 



q s \dE(q s )=0 



dEdQ 



Remark. The residue type formula for the third derivatives of F(T) implies that if we consider 
the reduced hierarchy with dQ = dp, then F(T) is a solution for the WDVV equation. 



Real-normalized Whitham hierarchies 

Real leaf and Whitham times We begin the study of the real-normalized case by in- 
troducing the following real-analytic coordinate system on A4 g (n,m): 



tl 



3? 



— iresp, z k QdE 
k 



tl = % 



— resp, z k QdE 
k 



k = 1, . . . , n + m, 
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and 



t x — V 



L h,a k — ° 



(IE 



dE 



b k 



QdE 



tX E,a k ~ $ 



l E,b k ~ ^ 



l h,b k ~ ^ 



dE 



dE 



i>k 



QdE 



k = l,...,g. 



tQ,a k ~ & 


I dQ 




/ dQ 




J a k 




■J a-k 


f Q,b k = K 


I dQ 




I dQ 


Jb k 


Jb k 



Note that all coordinates except t x ha , t\ b are just the real and imaginary parts of the 
corresponding complex analytic coordinates. A real leaf Vr(?t., m) C Ai g (n,m) is defined to 
be a zero set of the coordinates t v E ak , tE,b k JQ,a k ^Q,a k - Alternatively, Vr(ti, m) can be defined 
by the following Ag basis-independent equations, 



(16) 



dE 



= 0, 



dQ 



0, Vc G Hi(T g , Z). 



Therefore, Vr(ti, m) is well-defined as a real-analytic submanifold of both A4*(n,m) and 
M g (n, m). 

Before introducing Whitham coordinates on VR(n,m), we have to make the following 
important observation. In what follows we work with multivalued differentials that usually 
have the form fdcu, where / is an Abelian integral and the differential duj can again be 
multivalued. These differentials are well-defined only if we cut the surface T g . Thus, we 
always assume that some canonical basis B of T g is chosen, and we cut our surface along the 
representatives ak, bk of cycles in B. Sometimes these cuts are not enough and we need one 
extra cut 7. Let us introduce the following notation. 

Notation. By t g we denote the normal polygon obtained from T g by cutting along the a 
and b cycles. For each cycle a/ we denote by a+ its left and by a^ its right sides. Same 
notation applies to 6-cycles. The point of intersection of af and bf cycles is denoted by 
and the point of intersection of af and b g cycles is denoted by $ - By (r g ) we denote 

a Riemann surface obtained from T g by making a cut 7 from $ to Pi, and by (r s ) we 
denote its normal polygon. Note that topologically (r s ) ~ T g — P x . 



' cut 



The multivalued differentials that we consider are single valued on either f g or (t g ) 



cut 



and, as differentials on T g or (r fl ) , they can have jumps across the a and 6-cycles, and 
a cut 7. Note that for such differentials, the integral over a cycle depends not only on the 
homology class of the cycle, but also on the side of the cycle and the actual choice of a 
representative in the homology class. Thus, in order to make our construction independent 
of a choice of such representative in the homology class, for each differential fdu we introduce 
a corresponding cocycle [fdu] by the following procedure. First, we define [fdu] on basic 
cycles ak, bk by adding certain correction terms to § fdu, and then extend it to an arbitrary 
cycle by linearity. We use the following notation. 
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Notation. For any cycle c = a\a\ + • • ■ {3 g b g , by 

fdu = [fdu] (c) = ai[/du;](ai) + • • • + /3 fl [/dw](6 fl ) 



we denote the value of [fdu] on a cycle c. Note that [fdu](ak) and [/do;](&fc) still depend on 
a side of the cycle. We always choose a right side for a-cycles and left side for 6-cycles, 



fdu :-- 



/da;, <z> /da; := ® /da;. 

J[bk] M] 



However, it is convenient to indicate the side explicitly in the intermediate calculations. For 
consistency, we use same notation for usual Abelian differentials on T g . In addition, for any 
differential d/, we denote by / = f^' df the value of the corresponding Abelian integral / 
at the point $/ of intersection of the a~[ and fy + cycles. We also put 

n = » [/(<*>*)] , // = 9 [/($,)] . 

We now define Whitahm times on Vm.{ti, m) by 



9£ 



-1 

T 



and 



1 h,k — ^ 



T" 3 — c> 



res Pl z QdE 



QdE 

1 

dQ 



T k = $t 



-1 



•resp 1 z QdE 



k — 1, . . . , n + m 



[ofc] 



I h,k — ^ 



r r a c> 



Qd£ 



d£ 



[6*] 



rpa 

1 E,k 



rjib 

A Q,k ~ 



-9 



-9 



dQ 



[6fc] 



d£ 



[o/fe] 



Except for some relabeling, the main change occurs in the definition of T£ fc , T]£ fe : 

rpa A A rpx rpb it A rpx 

1 h,k ~ L h,a k L Q,a k rj kJ 1 h,k ~ L h,b k i Q,b k rj k- 

Prepotential The prepotential dS = QdE of the real leaf Vr(ti, m) is described by the 
following proposition. 

Proposition 1. The prepotential dS is a meromorphic differential on T g , holomorphic on 
T g — Pi, with pole of order n + m + 1 at Pi, and with jumps on a and b cycles. Near Pi, 



'm+n 



dS~d[J2(Tj+ y^lTj)t j + 0{z) ) + R s ^, 



where z — t 1 is our preferred local coordinate in the neighborhood of Pi, and 



(17) 



R s = res Pl dS = £ (T^T^ - 2*^) . 
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The jumps of dS across the a and b-cycles come from the jumps of Q and are given by 
(18) dS(P+ - P a -) = V^lT^dE, dS(P+ - P h ) = ^lT b Ejl dE. 

The corresponding cocycle [dS] is well-defined only on T g — Pi, and is given by 

AT b El E x I dS= I dS+^T a El Et 



dS 



dS 



The proof of this proposition is a direct calculation. 

Since dS has a residue at Pi, the corresponding integral S(P) = f dS is well-defined 
only on (f 9 ) cut - We choose an additive normalization constant in such a way that the regular 
part (5')_ vanishes at Pi. We can now consider cocycles corresponding to the differentials 
EdS and SdE, which are well-defined on T g — P ± . 

Proposition 2. The cocycles [EdS] and [SdE] are given by 

I EdS=l EdS + V^lT^St + Tq {T e i Ei - ^\T b E ^, 

J[a l ] Ja t Z 

I EdS=l EdS-V^T^ + T^T^Et + V^TT*,^, 
J[b+] Jbi * 

I SdE= I SdE-Mf^ + y/^lT^, 

J[a t ] Ja t ^ 

I SdE= I SdE - (io s ) b+ E l - V^1T E t 3-. 

Ml Jb t ' 2 

Moreover, the following "integration by parts " formulas hold: 





1 SdE 


= -3? 


/ 




Mi 




J[a-] 


3? 


(f SdE 


= -3? 


V 




Mi 




J Hi 



Whitham differentials Similarly to the complex-normalized case, Whitham differentials 
dVtA(E\T) = dT A dS(E\T) can be divided into four groups: 

• dQ x k and dVL\ are real-normalized meromorphic differentials with prescribed singularities 
at Pi, i.e., exactly the differentials of the real-normalized Whitham hierarchy, 

• dflfr k and dQ b h k form a canonical real basis in the space of real- normalized holomorphic 
differentials that is dual to our basis B, 

• dVL a E k and dfl b E k are meromorphic differentials with a simple pole at Pi and a dP-jump 
across a and 6-cycles, 
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• d£lq k and dVL b Q k are meromorphic differentials with a simple pole at Pi and a dQ-jump 
across a and 6-cycles. 

In the sequel we restrict to leaves V^(n, m) of a foliation defined by the level sets of Tg )fc , 
Tq k . Then all the differentials but d£lq k and dVt b Q k generate flows preserving the foliation. 
For this reason, we do not describe differentials dVLq k1 dVL b Q k in detail. We also need to 
consider corresponding Abelian differentials, which we always normalize by (Qa)(Pi) = 0. 

Proposition 3. For any Whitham time T A , the differential dVt A (E\T) = d TA dS(E\T) is 
holomorphic on Y g — P\. 

Proof. The only possible extra poles can appear when E(P) does not define a local coor- 
dinate, i.e., at the points q s such that dE(q s ) = 0. Assuming for simplicity that dE has a 
simple pole at q„ let E S (T) = E(q s (T)\T), and choose £(E\T) = ^/E - E S (T) to be a local 
coordinate near q s . Then we have 

dr A dS(E\T) = dr A (Q(((E\T)\T))d(e + E.) 

= ^=M^ 2 ^ = _ ( a r ^ s ) dQ(B | T ), 

which is holomorphic at q s . □ 
Main properties of the Whitham differentials are summarized in the following proposition. 
Proposition 4. 

• Differentials dVL\(E\T) and dfl k (E\T) are real-normalized, 

= VcefTi(r ff ,z), 

meromorphic on T g , with a single pole at Pi, where 

dQl~d(t k + 0(l)), dQ\ ~ d(</=U fc + O(l)) res Pl Q z k — res Pl Q k — 0. 

Differentials dQ% k (E\T) and dQ b h ^ k (E\T) are holomorphic on T g and form a canonical 
basis in the space of real-normalized holomorphic differentials dual to our homology 
basis B, i.e., 



3 


jdtll 


= 3 













3 



3 



/ dn a Kk 


— 8m, 


3 


/ dW Kk 
Jh 


= 0, 


/ dn b Kk 

-J at 


= 0, 


3 


<f d^M 
Jh 


= hi 
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Differentials dfl E k (E\T) and dQ b E k (E\T) are meromorphic on T g , holomorphic on T g — 
Pi, have a simple pole at P\ which is balanced by a single jump across one of the cycles: 



res Pl dtt a E k 



res Pl dVt E , 



1 



-1 



2tt 



T b Q , k = dE, dn E jp+ - P t 



a, 



■!SkidE(P c 



0.1 ) 5 



"k 



pz 



1 



rpa __ 



dE, 



dn b Eik (p+ - p- ai 
dti% k (p+ - p b -\ 



o, 

0, 



^18 kl dE{P bl ) 



differentials themselves are not real-normalized, but the corresponding cocycles are, 

0. 



3 



/ dsn,* 


= 3 




U[c] 




Ulc] J 



Duality and Riemann Bilinear Relations In the real- normalized case, the dual "inte- 
gral" operators § Ta are given by the following formulas: 



I dtt = $l [res Pl tdVt] , 



I dtt = M [res Pl V^TPdn] , 



dn = $i 
dfi = ^ 
dn = $t 



rpa 

E,k 



-1 

2^ 

-1 
2^ 

-1 
2^ 



dQ 



k+i 



dQ — ^R. 



1 

2^ 



dQ 



EdQ + GoTT a k (dQ) 



EdQ + CoiT b k (dQ) 



k+i 



where the correction terms Corr(<if2) are given by 

Conl(dS) = T% k Tq k + T hk Tq k , Con b k (dS) = —T^ k T b Q k + T^ fc Tg fc , 
Corr^OO =a TA Corr(c/5). 

We now have to establish relations flnp. 

Proposition 5. Whitham differentials satisfy the identity 



(19) 



(p dQ B = f dQj 
Jt a Jt b 



Proof. This proposition is proved by a direct calculation. We illustrate it by considering the 
following cases. 

First, let T4 = and T p = TJ. Then we have to prove the following identity: 
/ d£l) = K [res Pl Vdfflj] = 3? [res Pl t j dQf\ = I dfij. 
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But 



res Pl Pd&lj = - res Pl n)dt = - res Pl Q^dQ] - res Pl )_ 

= - res Pl VL)dVL\ - res Pl Vd(V)_ = - res Pl + res Pl (t j dQf). 

Taking real parts and observing that 



3? [res Pl tydn*] = 3? 



2tt, 



0, 



we obtain the desired identity. 
Let T A = T£ fe , T B = T« t . Then 



is just the usual Riemann bilinear identity for a canonical basis of real-normalized holomor- 
phic differentials and it is proved in the regular way: 



= K [res Pl Q a hk dfl a hl ] = 3? 



— ( - $ dttt k + <\> dn a Kl 



The most difficult identities to establish correspond to Ta = T E k , Tb = T E s and the 
like, since dVt a E k , dQ b E k have simple poles at Pi. In this case, let C £ be a small circle around 
Pi, $ e G C £ , and let 7e be a cut from $ to $ £ . We consider the integral VL a E k dVt a E s along 
the contour 



(20) 



First, we compute that 



J2(af + 6+ - a~ - b~) + 7 + - 7 ; + C £ . 



3? 



2tt 



1 9 



= 3? 



= 3? 



+ 



I tT Eik dE + 6 ks ElT b Qjk 

Ja s 

i- U Edn% s + £_ n% k dE + T^ k ((Q E J k + 5 ks El) 



— [ 4 Edn a Es + 



dQ% k ) El 



Note that 



(n a E J k + 5 ks Ei = (n E J + 9 



+ SksEl 



= (^,Jo = 9 fe($o )] = 3 fe(*o + )] 
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since the only non-trivial imaginary contribution for periods of dVt a Es comes from 6-periods, 
and in order for ft + -contribution not to be canceled by ^"-contribution, we need k — s, in 

'*0 



which case S 



On the cut 7 we have 



-6 ks El 



To evaluate J c Q E k dfl E s , we rewrite everything in terms of a local coordinate z in the 
neighborhood of P\\ 



+ (a + 0(z))dz, Ql k 



1 Q,k 
27T V / ^T 



Log(z) + 0(z). 



Then 



$£ ^Q,k^Q,f> L0f 



-4tt 2 



-d* + 0(e) 



J- n uJ- 1 



Q,k ± Q,s 

8tt 2 



(Log(*)) s 



Collecting all of the above together, we obtain 
= 3? 



3? 



r Q-^- a (47rv^lLog^($+)) + 4vr 2 ) + 0(e) 
; < v / /is i- 1 . ' Q,k_Qjl _|_ 0(e) 



z=i V a z+ 



^(^_ Edn% s + j _ n a E:k dE + o 



+ 

7, 



^E,kd^E,s 



3? 



1 

27 



n«dE 



as e — > 0, 



which proves our identity. All other cases are somewhat intermediate in difficulty to the 
cases considered above and are proved along the same lines. □ 

The third derivatives of F(T) are given by the following theorem. 

Theorem 1. 



d TATBTc F{T) = K 



q s \dE(q s )=0 



dVLAdQsdQc 
dEdQ 
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Proof. First we prove the following formula. For any two Whitham differentials dVt B , dVt c 
we have 



3? 



2tt 



\= I (d TA Q B )dQ c 
v— 1 ^af 



3? [res p 1 (5r A fis ) df2 



(21) 



+ 3? 



E 



res, 



q s \dE(q s )=0 



dVlAdVl B dVlc 
dEdQ 



To establish this identity, note that since dVtc is holomorphic outside of P±, the right hand 
side of our equation is a sum of residues at Pi and at poles of 8t a ^b- Using E as a coordinate, 
we see that 

d TA n B (E\T) = d TA d TB S(E\T), 

is holomorphic. So the only extra poles can appear at q s such that dE(q s ) = 0. Assuming 
that q s is a simple zero of dE and using £(P\T) = a/ E(P) — E S {T) as a local coordinate 
near q s , we have 



dQ 



n B (P\T) = n B (q s (T)) + ^£( q ,(T))t(P\T) + 0(e(P|T) 2 ) 
(3 tA )(P|T) = ^(q.{T)) -%£W + 0(1) 



2£(P|T) 



(^Q)(P|T) = ^fa(T)) + 0(1). 



At the same time, 



Therefore, 



dE(P) dP 



(2P 



P\T) + 0(1) 



and we have 



^fe(T))-^(P|T) + 0(l) 



res, s (a TA Q B )^ c = re S(?s ^M T ))^<^(P|T) 



= res„ 



dVtAdVL B dVt 



c 



(P\T). 



qs dEdQ 

The rest of the proof is a direct calculation, which we illustrate by three different cases. 
d TA {d Tl d T rF{T)) = U [res Pl ^d(d TA ^)] = 3? [res Pl ^d(d TA ^)] 



-3? [res Pl dr^dnf] = 3? 



E 



res 



q s \dE(q s )=0 



dvt A dn\dn) 

dEdQ 
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since #r A f^ is holomorphic at Pi and 3? 2 J^l $df. (^T A ^))d^l 

-1 



3? 



3? 



res. 



q s \dE(q s )=0 



dEdQ 



since resp t (dr A ^,i) dn M = °- 
d TA {d T%i d T%k F{T))=U 



-1 

2^ 



3? 



3? 



= 3? 



2tt 



/ Ed(d TA Q% k 
J K-] 

X ~r=i (dr A n% k )dn a Etl 

dQ A dQ a Ek dtt a El 



^z<f (dr A n E>k )dE 



res 

9s M£'(q s )=0 



since resp^c^f^^df^j. = an d (^t^e k) * s well-defined on f g . 



□ 



Reduction to Probenius Structures and Solutions to the WDVV equations In the 

real-normalized case we choose the marked variable to be T{. It correspond to the direction 
of unity, 



(22) 

In addition, we choose 
(23) 



d T x = e ~ dQ[. 



dQ = V-ldQl 



and truncate the hierarchy by considering T/, T\ only for i < n. Then we have the following 
theorem. 

Theorem 2. Given the above assumptions, 



(24) 



dT{T A T B F(T) - < 



ij 
n ' 


dQ A = dQj, dflB = 


dty, and i + j = n 


ij 
n ' 


dVl A = dVL\, dVL B = 


dQ v j, and i + j = n 


1 

2tt' 


dVt A = dVL a E k , dVL B 


= dK,k 


1 

2tt' 


dfl A = dQ b E k , dVLs 







otherwise. 





Therefore, the coordinates are flat and F(T) is a solution to the WDVV equation. 
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Proof. Let us rewrite the formula for the third derivative with dQ = \/—ld£L\: 



dT{T A T B F{T) = 3? 



3? 



res„ 



q 3 \dE(q a )=0 



dVtAd^B 



dVL^dVlB dQsAd^B 
-res Pl 



2tiV-1 /af g V-l^ V-ldE 
First let us see what is happening near the puncture P\. Since dE ~ d(z~ n + 0(1)) near 



Pi, ^| ~ has a zero of order n + 1 at P\. It is clear that for the truncated hierarchy 
the only differentials that can contribute to the resp 1 (*)-term are d^l\ and d£l\. Assume that 
near Pi, 



Then 



dVL A ~ d(sAZ-' 1 + 0{l)), dQ B ~ d(e B z- j + 0(1)), e = 1 or V^l. 

dft A dft B = (^(-i)*"*" 1 + 0{l)){s B {-3)z-i- 1 + 0(1))(^) 2 

= {e a e B (tj)z-^- 2 + O^" 1 ) + O^ 1 ) + 0(l)}(dzf 



3? 



res Pl 



= TdE 



3? 



V — l^A^B 

— n 



A — —J\ A 



Now let us consider the boundary term - — = <f Q p ^j=f^§ ■ For it to be non-zero, at 
least one of dflA,dil B should be dVL a Ek or dVL b Ek . For example, let dQ A = dVL a Ek . Then, if 

dn B ^dn E>s ,dn b EjS , 



dE { ai ai) 



-IS, 



dQ,% 



E,k 



kl 



dE 



p: 



o 



dVL a Ek dVL B 

27TV^T Jgt v y^ldE 



and if dfl B = dVL a E s or dVL B = dQ b E s , 



2tt^T 



I, dn B = dni k 

otherwise 



3? 



2tt 



d^E,k^E,s 

IdE 





2ttv^T f a t g 


= 3? 







dn% k dn b 



IdE 



0. 



All the remaining cases are similar to the ones considered, and this concludes the proof of 
the theorem. □ 
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